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The hygrothermal transfer is very important for the design of a building 
envelope for thermal comfort, economic and energy analysis of the building 
envelope. The lack of reference materials on models of moisture and 
temperature behavior in the building, including wooden walls, is a challenge. 
This paper reviewed the hygrothermal transfer models for building walls. 
Energy and mass conservation equations with boundary and input conditions 
were presented in this paper for concrete, bricks, and wooden walls. The review 
showed the presence of mainly physical-based models, while there is a dearth 
of data-based models. The influence of the type of wall, orientation, thickness, 
the density of the material, and climatic variations on the temperature and 
moisture evolutions within the building materials influenced the model 
mechanisms. Future research gaps should include shrinkage influence on 
hygroscopic materials like wood due to their behavior under ambient 
conditions. Data-based models should be explored too. 


1. Introduction 


Heat and moisture transport has 
simultaneously in building envelope. The relative humidity of 
the air indoors can affect the micro-climate of the building 
envelope and, by extension, the energy consumption. While 
there is aclamor for energy-efficient buildings [1], this should 
be achieved in consonance with the material that has high 
moisture buffering capacity to avoid moisture damage within 
the building. Some researchers have estimated that building 
alone consumes about 36 -70 % of global energy, generating 
close to 50 % of global greenhouse gas emissions [1-3]. This 
energy is expended right from the operational phase of a 
building, material extraction, production, 
transportation, and the end of the life span of the building [1- 
9]. Therefore emphasis now is to cut down energy utilization 


been 


construction, 


in a building to reduce greenhouse generation on the 
studied environment, acidification of the environment, depletion of 
the ozone layer, global warming, abiotic resources reduction, 
eutrophication, etc [1,10,11]. Countries are adopting a 
different green strategy in building, hoping to cut down 
energy consumption by 20- 42 % by 2050 with a 35 % 
reduction in greenhouse emissions [12, 13]. While some have 
considered the entire structural envelope, others have looked 
at the materials for construction or various components of the 
structures ranging from the walls, the floor, or the roof 
envelope [9, 14, 15]. The adoption of hygroscopic materials 
throws up the issue of moisture adsorption for these 
hygroscopic materials. Condensations inside the building 
envelope are most common due to variations in temperature 
and humidity indoors and outdoors the building. The 
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presence of moisture within the building envelope will cause 
thermal discomfort and can corrode the metallic structures 
within the building and also make the indoors moldy [16]. 
This can lead to structural degradation and failure. Therefore 
the study of heat and mass transfer between the porous 
materials interface is done to elucidate their thermal 
performance and strength. Understanding the physics of this 
physical process is important to be able to predict them. The 
discontinuous moisture profile of two porous materials at the 
interface, because of their hygroscopic characteristics, has 
been used to predict the temperature and moisture gradients. 
Therefore, it is obvious that the nature of the materials used 
in building construction affects the temperature and moisture 
variations within the building as a function of the ambient 
weather changes and, by extension, the energy required for 
cooling or heating satisfaction [9, 17]. This will determine the 
magnitude of humidity, heating, and cooling required [18, 19]. 

Materials like wood, metal, steel, concrete, pozzolan, 
glass, bricks, and other multi-layered composites have been 
adopted for building walls [20-25]. Osayintola et al. [26] and 
Lelievre et al. [17] classified building materials as classical 
and hygroscopic building materials. Lelievre et al. [17] 
further sub-classified the building materials as bio-based and 
non-bio-based materials. The authors stated that bio-based 
materials like hemp concrete have potential low carbon 
emission, good thermo-hydric properties, and moisture 
buffering properties [27-30]. In several cases, walls can be 
made of more than one material layer, considering building 
insulations and cement plasters in some walls. Insulated walls 
(Figure 1) will regulate the heat and moisture transfer within 
the building envelope. Therefore, experimental and 
numerical studies have been carried out to study the physics 
of heat and moisture transfer of various materials in response 
to variation in weather parameters, which includes solar 
radiation, temperature, moisture, and relative humidity. 
These models are developed and present as single or 
multidimensional cases in the literature [31-47]. These 
models are resolved using finite elements, finite control 
volume etc, in a steady and non-steady numerical scheme. 
They are majorly predictive models to determine the 
temperature, relative humidity, and moisture condition of the 
inside of the building as a variation to the ambient conditions 
[32-34, 36-43]. The various studies used the thermo-physical 
properties of these materials, the nature of airflow, the 
dimension of moisture transfer, and the meteorological data 
of each area to develop the simulation codes to make their 
predictions with good result. The vapor adsorption and 
desorption isotherm is controlled by the main adsorption 
isotherm with the occurrence of hysteresis in the sorption 
curves [17]. Although the nature of hysteresis is yet well 
established, some authors included it in their modeling 
approach, while others neglected it. Simo-Tagne et al. [9] 
stated that studies like this would help in the selection of 
materials to satisfy the desired heat load at minimum 
dissipation of energy with less environmental impact. Some 
researchers have tried to review these models, though they 
presented only computer-based model tools, like UMIDUS, 
Wuf, MATCH, DELPHIN etc [46, 47]. Prior to that, the 
Canadian mortgage and housing cooperation in 2003 
reviewed and showed that about forty five hygrothermal 
transfer models were in existence, which increased to about 
57 models by 2008 [47]. The use of a dynamic coupled co- 
simulation approach for forecasting the hygrothermal 
behavior of building envelopes was discussed by Ferroukhi et 
al. [46]. However, with the development of new material and 
the creation of different kinds of boundary conditions and 
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different interactions with the environment considered, new 
numerical models on hygrothermal transfer are presented 
and validated with experimental data. Busser et al. [78] towed 
this line in presenting the recent trends in the experimental 
validation of hygrothermal transfer models, but the specific 
equations were lacking in their presentations. 

Prior knowledge of energy behavior is important to 
design a building for optimum comfort, taking into 
consideration all the heat loads. Moisture and temperature 
changes in the building are coupled together to study this 
moisture and temperature gradient. This requires the 
development of modeling tools for optimal condition 
prediction for different kinds of walls. Due to the effect of 
greenhouse gas emissions generated from the production of 
non-bio-based walls, interest is shifting to environmentally 
friendly walls. Wooden walls are now of interest, especially in 
Africa, when the cost of non-biobased walls is also considered. 
This review is an updated review of hygrothermal transfers 
with the addition of research work on wooden walls lacking 
in other previous reviews. Therefore these will bring up to 
date various reviews conducted with the same theme. The 
review will single out each model and discuss its pros and 
cons. 


2. Methodology 

This review involves searching and requesting available 
open literature on hygrothermal behavior for different walls. 
Emphasis was on the literature on models that have not been 
reviewed, although where the models are based on existing 
models, the old models were presented, and the new models 
can be discussed under them. Subsequently, the models were 
separated into those developed and validated with concrete 
and bricks and those validated with wood. After reading and 
deducing the mechanisms governing the behavior of the 
hygrothermal transfer and important equations and results, 
this is presented and discussed. Figure 1 gives the flow chart 
of the review methodology. 


Energy and moisture transfer in buildings: A review of 


hygrothermal transfer models through the porous 


Yy 
| Literature sourcing ) 


Models for concrete Models/tói woods 


and brick walls 


Discussion 


Figure 1. Schematics of the review process 
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3. Models for concrete and brick walls 

Several modeling and simulation investigations exist 
with a different approach for concrete and bricks based walls 
in literature. To reduce energy consumption or improve 
hygrothermal transfer in buildings, researchers have 
developed different kinds of wall composite. Walls can be 
single, double or multilayer. A multilayer wall is a wall made 
up of more than one porous material. Concrete (bio-based and 
non-bio-based), bricks, cement plasters etc can be used to 
develop a multilayer wall. These walls have also been 
numerically simulated to obtain different varying conditions 
for indoor and outdoor conditions, and the results of the 
models were validated experimentally. The various models 
are as follows. 


3.1 Xingguo models 

Xingguo et al. [20] modeled a multi-layered porous wall 
made of cement mortar, red bricks, and cement plaster in the 
southern Chinese city of Hunan. The model developed was a 
modified one-dimensional transient hygrothermal model 
with temperature and humidity as the driving potential. Two 
key transport equations (1 and 3) for mass and heat 
transport, respectively, were solved using the finite element 
with established boundary’ conditions and the 
thermophysical properties of these materials. The maximum 
temperature and humidity difference obtained by the 
researchers were 1.87 °K and 11.4 % for indoor and outdoor 
conditions. 


aw _ Ws 


W , OW, aT 
at Pm 


Pm a? 
(DyRyTmpa T W; Dw) ax? t aT at (1) 


The resultant boundary conditions for the above equation 
were given as follows. 


ow 
—DyRyTmPa oe Am (Woo = Wurf) (2) 


For heat transport, the following governing equation was 
used. 


oT 0°T a*w 
(PmCym) a Ka + hfgDvRyTmPa Ox (3) 


The resultant boundary conditions for the above equation 
were given as follows. 


OT 
-K z hm(To Ta Tsurf) + Qraa F hfgms (4) 


3.2 Lelievre model 

Lelievre et al. [17] combined two sub-models of 
Pederson and a phenomenological model from Mualem II to 
develop a numerical simulation model for multilayer hemp 
concrete. Hemp concrete is usually coated with plasters of 
different levels of permeability inside and outside, and the 
thickness is non-homogenous. The model developed, which 
accounted for phase change and hysteresis, depended on the 
temperature and moisture transfer as a function of the 
hygrothermal properties of the hemp. They gave the energy 
and moisture conservation equations as follows. 


ps(Cp,s + WCp1) = = -V x (—AVT) + -V x (Df Vo + 
DEVT) x (ly + (Cps — Cpa)(T — Tres) ) (5) 


ps0 £ = -V x (-(D? + D?YT)Vo — DEVT) (6) 


The sorption capacity of equation 6 was deduced using two 
hysteresis models in sorption and desorption phases from 
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Pederson (equations 7 and 8) and Mualem (equations 9 and 
10) as follows. 
= B(w-Waa)4Odes+(W—-Waes)*Oaa 


(Waes—Waa)4 


(7) 


Oad, hys 


= (W—Waa)4 Oaes+ C(W—Waes)“Oaa 


(Waes—Waa)4 


(8) 


Pada nys 


Using Mualem models, the following sorption equations were 
deduced. 


Waes,nys(P) = Wj es (ws — Waa (y)) (Waa (9;) — Waa (9)) 
(9) 
(Waa (9;) — Waa ()) 
(10) 


Wj-Wi 


(Waa(¥j)—Waa ()) 


Waa,nys(P) = Wj —- 


Validation of the above models showed that using sorption 
isotherm from Mualem gave a good agreement between the 
experimental and predicted results, while the model of 
Pederson was off the mark. 


3.3 Djongyang model 

Djongyang et al. [40] presented a hygrothermal transfer 
model for porous building components, which they validated 
with earth bricks wall. They considered a plane geometrical 
shape and the influence of inter-tropical conditions with 
variations in latitude for the three cities of Cameroun. In 
solving the numerical equations, they considered the periodic 
solution approach and validated their model with two works 
of Menghao et al. [48, 49]. The model developed was one 
dimensional in which liquid water and air and water vapor as 
a single binary gas mixture was considered. The conservation 
equations for mass and heat transport were taken from the 
equation of Luikov, which has been used by other researchers 
[50] as follows. 


du(x,t) 0? u(x,t) 0?T (x,t) 

ax Om jy + amô = ae (11) 
aT(x,t) _ a? T(x,t) du(x,t) 
ae ae + ef a’ 0<x<l (12) 


Te boundary conditions for the two equations above were 
defined as follows. 


OT (x,t) 
—Kgout = F hout (Tx=0 = Tout) + Aout = 
Ox x=0 
Eout) Ux=0 = Uout) =0 (13) 
u(x,t) OT (x,t) 
Kmout = i + Kmout Sout a a + Aout (Ux=0 = 
Uout) = 0 (14) 
OT (x,t) 
—Kgin —— =o + hin Tx=0 = Tin) t Am > Ein) (Ux=0 = 
Uim) = 0 (15) 
u(x,t) OT (x,t) 
Kmin = + Kmin Sin x + Qin Ux=0 = Uin) =0 
x=0 x x=0 


(16) 


Equations 13 and 15 represent the heat balance with the 
three components of the quantity of heat exchanged for 
outdoor and indoor, the convective heat transfer and 
evaporative flux, while equations 14 and 16 represent the 
mass balance with the components of moisture gradient, 
temperature gradient and the convective flux exchanged 
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between the ambient and the surface of the materials. 
Djongyang et al. [40], considering the length of the day and 
the declination of the sun, converted equations 11 and 12 to 
oscillatory (periodic) equations using the Fourier series 
method as presented in equations 17 and 18, which they used 
in their validation adopting a periodic approach 


07un(x) inw 3? Tn (x) 

3x2 — “Gm e" +6 3x2 =0 (17) 
0°T (x) inw inw 

Ax? cara in +ep | Un =0 (18) 


Using a numerical approach, they solved equations 17 and 18 
to generate the temperature and moisture profile of the 
material under variable external conditions. Validation of the 
results using thermo-physical properties of earth bricks gave 
good results. They showed that latitudes affect hygrothermal 
transfer. 

Menghao et al. [49] also used the same one-dimensional 
Lukoiv equation presented in equations 11 and 12 in the 
modeling of a hygrothermal transfer for a fibrous slab, but 
they added the effect of adsorption and desorption heat, 
which is also one of the source or sink terms in coupled heat 
and mass transfer equations. They assumed a localized 
thermodynamic equilibrium between the fluid and the porous 
matrix in presenting the equations as follows. 


oT 3T a 

Cops, = kaa + Cmp (eh + 1) 5, (19) 
ð a? 3T 

CmP M = Dm oer +Dmô = (20) 


The latent heat of vaporization will be integrated as part of 
the energy balance at the building materials boundary, which 
will be affected by mass diffusion due to deviations in 
temperature and moisture content [40, 51]. The boundary 
equations were given as follows. 


Tee -o Z (T (0,0) = T1) + fihi — €)@n(0, t) — m) 
(21) 
-KZO = aT (Lt) — Ta) + Bahy (1 — E)M(L t) — 
mi (22) 
Din OT (x,t) a kmô aT (x,t) = By (m(0, t) a mı) 
ðx Ix=0 =0 
(23) 
-Dm OT (x,t) — km 6 — aT (x,t) = Bo (m(L, t) = mz) (24) 
ðx İx=0 ax =0 


In the resolution of the equations using Laplace 
transformations, dimensionless terms were introduced to 
equations 19-24 to change it to dimensionless form. 
Validation of the equations was done using the experimental 
setup in Figure 2. The porous material is a multilayer material 
whereby that one side is permeable and exposed to outdoor 
humid conditions while the other side of the wall, which is 
impermeable, is subjected to cold temperatures. Validations 
of the numerical simulation model, although showed the same 
trend as the experimental data looking at Figure 3 plotted by 
the authors it shows that there is still temperature variations 
which implies that there is still some omitted parameters like 
hysteresis, which is yet to be understood by researchers 
because of non-homogeneous nature of some porous 
materials; however, moisture profile was more accurately 
predicted. 
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Figure 2. Experimental set of the testing rig for model 
evaluation [40, 51] 
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Figure 3. Comparison of the numerically simulated 
temperature profiles with the experimental result after quasi- 
steady state [40] 


3.4 Simo-Tagne model 

Simo-Tagne et al. [9] presented a numerical model for 
predicting the hygrothermal transfer for concrete walls for 
outdoor conditions in sub-Saharan Africa. The model took 
into account all kinds of bound and integrated the type of flow 
in the boundary conditions, which is not common in other 
established models. They also stated that all water type 
occurring in the material is modified during the moisture 
transfer; therefore, they presented the mass conservation 
equations for liquid water, bound water, and vapor phases in 
equation 25-27, respectively, as follows. 


ae +77, = (25) 
að X Cad = 
OOD A T Jas = —Kas (26) 


The contribution of water vapor to the mass balance 
equations is given by: 


Heo) 6 7. (pg ¥y + Ig) = Kas + Ki (27) 


Pg; is the flux characteristic of the movement of the vapour 


phase. The solution of equation 27 gave equation 28 as 
follows. 


> 


fa] 
20S) +7 (pgVy +J +Jas +Jg) =0 (28) 
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The heat energy balance was given by equation 29 as follows. 
aT, eeg aT, pe 

[ec Z+ Vir| + [pc +0. 7r] + [pc = +7. 7r] Ee 

KıL + Kgs(L + Ey) = 0 (29) 


The boundary conditions are given in equations 30 and 31 as 
follows. 


aT 
(A+ oPs(Ep + L)Dy) z 
x=0 


OH 
[-ps(Ep + L)Dy = — 


heere Tent = T(O, t)) F PwLhm,ext(Xeq(Taext HRaext) m 

H(0, t)) +Ggioi (30) 
H OT 

[-D. a &tDy || = Henexe(Xeq (Taexts HR ext) B 

H(0,t)) (31) 


The above equations were discretized with finite differences 
using the Crank-Nicolson scheme, and the solution is with the 
Gauss-Seidel relaxation iteration method using Fortran 90 
language. Validation of the model was with hemp concrete 
with a low mean relative error. 


3.5 Philip and De-Vries model 

Philip and De-Vries proposed a model applied to study 
the heat and mass transfer through a porous medium using a 
building wall based on the thermodynamics of irreversible 
processes [78]. This model has been the bases for some other 
models, which will also be discussed in this section. This 
model took into consideration the effect of temperature and 
humidity gradients, and also, the altitude of the location was 
taken into account. The medium is assumed to be stable and 
homogenous. Kevin’s law was applied at every point of the 
medium, and the hysteresis effect between adsorption and 
desorption phenomena is neglected. This model is given as 
follows. 


Zyta Stet = div (D ggradx + Drgradt) “he 


ply i = div (agradr) + p,Ldiv (Da „gradX + DrygradT -= a% - c2) + 
pG (D gradX + DrigradT — K) gradT + ply (D qvgradX + DrogradT) gradT 
(32) 
kıg Pe ay 
ere a pw Cae (33) 
= oy Pus (OY Mygy 
Dr, =K (2) .p = fD “i 7 w exp (=2 ), 
Dq = Dq + Dav (34) 
P  (My\* LPys Mygy 
i P-P, a a exp ( RT ) r 
P My 2 IP us ( (OY Y 
Pp, = Pı (Ga. 3 HR (35) 
U-qvDqv PyMy 
T tit Dry, a a pik? 
U-qPry _ U-)Py 
= -7 36 
fda pıRT (36) 


Neglecting the heat transfer due to the mass transfer and the 
local temporal variation of the condensed water content in 
the vapor state, the simplified model of Philip and De-Vries is 
obtained, and it is well described by the experimental data 
obtained by Larbi in 1990 [78]. 
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Ox 
ot 


= div (D; gradxX + Dr,gradT) = x 


pC, me a = div (agraar) + pıLdiv (D qa gradX + DrygradT ) 


(37) 
The model of Umidus [79] assumed that the humidity is 
transferred through the wall only in vapor and liquid forms. 
The liquid form is transferred by capillary, and the vapor form 
is transferred due to the partial pressure of the vapor. Thus, 
the model in one dimension is as follows: 


ae =o (Pr gn) * ax (Di Se) 


a ar 

2 (Da 2) + 
Pin ôT _ ð (,0T al rv 55) 
PolCpo + Cy 2X) T= 2 (AT) + Ly 1 


ð Ox 
az (Dw 3x) 
(38) 
The boundary equations are given by: 
Humidity: 
ôT Ox 
—p\ (Dr ax + Dy ax) y=08 a en eee u Prese) (39) 
—Pı (Dr ax T4 Dy Z =Li = hy i(Pveai i= Poeni) (40) 
Temperature: 
ôT OT Ox. 
= = Ly pi (Dry ax + Duy or = = Are(Tae = Tia) + 
Lyhm,e(Pve,ae = Piese) (41) 
-4% — Lypı (Dr — +D (= hri(Tsi — Tai) + 
ax Lyi ( tty ae ToN TiN" si Tai 
Lyhui(Pveai — Pve,ai) (42) 


Using this model to study the heat mass transfer through the 
wall building with and without a perfect contact on a double- 
layer wall, Le [80] shows that the assumption with a perfect 
contact can produce some errors during the determination of 
the physical parameters influenced by the humidity. When 
the contact is real (not perfect), the model shows that the 
energy consumption is reduced by 10%. Using the 
experimental data taken from the literature, Le [80] validated 
this model. Using the same assumptions by Philip and De- 
Vries, the model of Duforestel is given by Simo-Tagne [78] as 
follows. 


Op, arp. OT 
ar ar wre (Üm +L) at 


v4 p PLKipoT Deco 
c= aiv(( Kp, + A wear) gradpy pikt grar ) 


P 
L(>+ 
L)\ op, | 1 aTPv i or : G 
oa at (e H gore hm im : :) at æf ‘ ee 
+Agra 
(43) 


Ky) gradp, 


(-hmar + hm 


This model is easy to use than the one of Philip and De-Vries 
and can be applied to non-hygroscopic material of 
construction, but the fact it did not take into account the 
gradient of moisture content is often the reason for the 
difference between experimental data and numerical data. 
Based on the same assumptions of the model of Philip and De- 
Vries, Luikov presented his model with other parameters. 
This model gives satisfaction when the partial pressure of the 
gas phases is uniform during the heat mass transfer. Luikov’s 
model is given by Simo-Tagne [78]. 
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= div [am (graax + égradr)| 
(44) 


ply Sr = 


= div (agraar) + eLŽ 
The difficulty faced in using this model is the difficulty in 
obtaining the conversion factor. Though it is mostly obtained 
using the inverse method, which is cumbersome. When the 
wall is constructed by the mortar slab, Saadani et al. [81] 
show that Luikov and Philip and De-Vries models can give 
similar values of temperature evolutions with some 
differences in the evolutions of moisture content, using the 
time steps between 0.1 and 1 s. Nitcheu et al. [82] used 
Luikov’s model with satisfaction to estimate the numerical 
heat mass transfer through the wall building constructed in 
earth bricks stabilized with thatch fibers. They used the finite 
differences method and implicit scheme of Crank-Nicolson to 
generate the numerical results, but the thermophysical 
parameters of the medium were considered constant during 
the process. 


3.6 Whitaker model 

In 1977, Whitaker presented a model based on the 
morphology and the average thermophysical values of the 
building wall. The mass transfer of water in gas and liquid 
phases were separately written. This model is very difficult to 
use, but Thaï used it with satisfaction in 2006 [78]. This model 
is given by: 


= 
ð Ew 
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Using his model, Nicolas studied the influence of the variation 
of the total pressure of the gas on the phenomenon of 
imbibition of cement. He obtained that the variation of the 
moisture content is influenced by the variation of the value of 
the total pressure of gas. The high pressure of the gas 
decreases the rate of transfer of humidity in the medium. This 
model permits the study of the effects of the hysteresis 
phenomenon during the heat mass transfer through the 
humid building wall [78]. 


4. Models for wooden walls 

The impact of wood as a building material for walls has 
been explored greatly in literature. An investigation has been 
carried out on the effects of wood on thermal comfort in terms 
of relative humidity, temperature, background noise levels, 
and CO2 concentrations [53-59]. Therefore, wood is energy 
efficient with lower CO2 concentrations. Wood, as a 
hygroscopic material, absorbs moisture from the 
environment, which has a direct and indirect impact on room 
conditions and thermal comfort [60]. Therefore, the type of 
wood, thickness, moisture isotherm etc. has been factored in 
the predictive models to describe the hygrothermal behavior 
of wood for indoor and outdoor conditions. 


Pw ae + Pw < Vy Hone 


2 (eg < py >9) + V(< Pa >9< V > 


ð 
24505 = = <T> 


In 1992, Nicolas proposed a model that took into account the 
capillary and hysteresis effects with the temperature gradient 
as a motor term. The Kevin’s law has been modified, and the 
pressures of water in saturation state, respectively in porous 
medium and in the free medium, have been differently 
written. The model obtained is given by: 


pCy (= + V. grad’) = div ig 
—((-¢,)T-L o) = -tn (2) 


Metu + div (bapava) =0 (46) 
a) Sie (pvb Vp) + + epa + div (p.biv1) = 0 
acre 4 div (bev) == = ae” n(2*) 

With : 


baVa = — (K; + Pa Kay) gradpa — (Kg — Kqy)gradp, (47) 


byVy = —(K, - Kav )gradpa — (K; + Pa Kay) gradp, (48) 


by, = —K,grad (pı + PT h(b1) + pig) (49) 
aRT RT dh(b,) 
Pa = PES py = 2 pi = Pa + Po + Poy (50) 


-7 |< po >9 Desp- V r (2)] 


<Pg>9 


(45) 


+ | Pwlew <Vy > +< pg >9< Cp >9< Y, >|V < T > +Ahy < M >= Pirar) 


4.1 Luikov models 

Due to the flexibility of Luikov models, it has formed the 
basis for modeling the coupled heat and mass transfer for 
porous material independent of its hygroscopic nature. The 
model accounts for all the bonding water simplistically 
without restricting the water transfer mechanism [61]. 
Therefore, these models can be applied for one, two, or three- 
dimensional coupled heats and mass transfer in woods. 
Younis et al. [61] used this model to predict the temperature 
and moisture transfer in a wooden slab. After defining the 
entire dimensionless variables, they used the finite element 
method, which operates in MATLAB, to solve the partial 
differential equations in three dimensions. They concluded 
that different dimensionless numbers (Luikov, Kossovitch, 
Posnov, and Biot numbers) in the coupled heat and mass 
transfer equations affected the overall heat and mass transfer 
behavior of the wooden slab. Comparison of this model with 
analytical and experimental data showed closer agreement 
with the analytical model rather than experimental data. 


4.2 Osayintola model 

Osayintola et al. [26] used a combined Knudsen and 
Fickian diffusion and neglected thermal diffusion in the 
modeling of heat and mass transfer evolution in spruce wood. 
The energy and mass conservation equations were defined in 
equations 51-54 as follows: 


0 +h =0 51 
p (51) 
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A(Pvég) _- _ 9 apy 
ot TE Ox (Des ee) (52) 
ar... _ a apy 
PCpeff 5, + Mpg = 5, (kerr oe) (53) 
ta = —p0% (54) 


The above equations were discretely solved, and a stable 
solution was provided with a relaxed, Gauss-Seidel iteration 
method. Discretization was by finite difference method with 
2d order accuracy for the implicit scheme. Although the 
authors suggested that the presented diffusion model can be 
used to compare the experimental results and set 
benchmarks for similar materials, they did not show this with 
their experimental data. Rather they only fitted the sorption 
data with moisture content, thermal conductivity, and water 
vapor permeability with good results. 


4.3 Simonson model 

Simonson et al. [69] presented a numerical simulation 
model to predict the moisture transfer (indoor climate) 
between the indoor air and the structural materials for 
wooden buildings in Belgium, Germany, Finland, and Italy. 
The model was developed to predict the temperature and 
relative humidity and applied to measure the comfort 
thresholds of the occupants in the building. The analysis of the 
model showed that the permeable interior layer is more 
satisfactory than when the interior layer is made with a 
water-resistant layer. The model is presented following the 
Ficks law of diffusion as follows: 


qM = —kd(u,T)Vpy — PoDw(u,T)Vu + Vapy + Kpwg (55) 


Equation 55 is adopted from IEA ECBCS Annex 24 ‘Heat, Air 
and Moisture Transfer in Insulated Envelope Parts presented 
in detail in Hens [62]. The model took into consideration all 
the energy components of the moisture transfer process for 
adsorption, desorption, condensation, evaporation, freezing, 
and thawing. The conservation equations were solved 
simultaneously to predict the variable indoor conditions 
under different experimental data obtained in the field [64- 
69]. 


4.4 Talukdar model 

Talukdar et al. [70] used a numerical model to set a 
benchmark for 1-D transient heat and moisture transfer 
models of spruce wood as a building material. Although the 
partial differential equations presented in equations 25-28 
was used to set the governing energy conservation equations, 
they used an analytical approach to set the moisture 
penetration depth through the woods using equation 56 and 
57. For a semi-porous material, the analytical equation 
developed for vapor transport is given in equation 30 as 
follows. 


(Pv-Pvi) _ ( x ) _ | (= k | x 
(Poo-Pvi) erfe Vamefft Deff . Derr 

x hmy&mefft 
|erfe ( arer + Der; ] (56) 


(Pv5m—Pvi) = 


The obtained results show that increasing the air velocity 
increases the temperature, relative humidity, and moisture 
accumulation within the plywood. 
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4.5 Watt model 

Watt et al. [71] proposed a two-dimension model for a 
light timber wall with an air barrier with changes in air 
tightness in a Swedish. The model was used to predict the 
moisture accumulation and mold growth within the building 
envelop. The magnitude of moisture accumulation is higher 
behind the outdoor air-tight layer of the simulated wall with 
the interior of the wall unsealed when compared to the sealed 
inside. 


4.6 Simo-Tagne model 

Simo-Tagne et al. [9] numerical model presented in 
equations 25-31 above was also used to predict the 
hygrothermal transfer for different wooden walls for outdoor 
conditions in sub-Saharan Africa. In this case, the model took 
into account the bound water presents in the wood with the 
integration of the flow pattern in the boundary conditions. 
Validation of the model was with Norway Spruce wood with a 
low mean relative error. Further analysis of the model 
showed that less dense wood with increased thickness 
provided better thermal comfort. Therefore, the nature of the 
wood and climatic factors is an important consideration that 
affects the hygrothermal transfer in woods. Simo-Tagne et al. 
[72] also presented a numerical simulation model for building 
walls in Nancy, France. The same approach and solutions in 
Simo-Tagne et al. [72] were adopted, but they integrated the 
Dufour and Soret effect. The driving potential for the heat and 
mass transfer was temperature and moisture gradient. The 
model showed the negligible influence of the wooden 
structure, while thickness is important in canceling the effect 
of ambient conditions. Simo-Tagne et al. [74, 75] present a 
novel model of heat mass transfer through wooden material. 
The model-based initially on the description of each type of 
water (bound water, vapor, and free water); equations 
obtained are given by: 


aw a aw aT 
at ax (Dun ae Dyr S) (58) 
ar_ a aw aT 
pl ~~ = 2 (Dr + A+ Drr) =) 
With: 
= p, — 2E (5x) 2Pe q PaPa 2Co 
Dny = Dy Ps a aw ` ps(1-Cg) aw (59) 
= D, — RE ie) 2e q ton Oe 
Dur = Dr Ds E) ar ps(1-Cg) OT (60) 
= Breas Oly _ Kr) OPe 
Drr = 1-Cy OT Epik ($), aT (61) 
Dry = See oe — Epik (1) 2 62 
TH 1-Cg êw tm Vu), aw (62) 
EpHR ôXe 
Tee a aHR " (63) 
The boundary equations are given by: 
If x=0, thus: 
aw oT 
“= 0 ao 0 (64) 
If x=+e/2, thus: 
aw 
—Dun am hm(W — Xe) (65) 
aT aw 
-(4 F Drr) ax = h(T = Tair) + PsLDuy Ox. (66) 
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Using the real variations in thermophysical parameters 
presented in equations 47-55, Simo-Tagne et al. [75] showed 
that the model is flexible and will allow taking into account all 
properties of the wood types and the movement of all types of 
water (bound water, free water, vapor of water). However, 
they suggested that to adopt the model, the hydric 
diffusivities have to be experimentally determined for each 
wood. Generally, to have an accurate model for all the heat 
and mass transfer model equations, the sorption and 
desorption isotherm for all the materials have to be 
accurately determined, and also the accurate boundary 
condition is defined [75-77]. These values and equations 
mainly differentiated most of the presented equations in this 
review from each other. 


4.7 Berger model 

Berger et al. [73] proposed a numerical model using 
Scharfetter-Gummel scheme combined with a two-step 
Runge-Kutta approach. Three phases were distinguished: 
water vapor, liquid water, and dry air. The moisture mass 
balance was given as follows: 


a 

= (Wy +w) = —D. Cx + Jot) 

: at (67) 
z Wr + Waa) = V. ov + Jeda) + lev 

The energy balance was given as follows: 


ôT 
(CoPo + CyWy + Ciwı + CaaWaa) a —ViJg- Tyley — 
V. (CT Yew — V-(CT Je — V. CaaT Ve,aa (68) 


Where the volumetric vapor source levis given by: 


= +P Jey (69) 


l 
cy R,T? ôt 


The heat flux was expressed as: 


T 


Ja = —À VT + (Cowy + CaaWaa) a-o) 


(70) 


Vis the vapor velocity taken equal to the air velocity and given 
by: 
kvda 


The flux of water vapor was given by: 


Py 
R,T 


Jev = —k, VP, + -5V (72) 


The flux of dry air was given by: 


J _ WvtWaa 
cda ~ (1-0) 


V = Jew (73) 


The flux of liquid water was given by: 


Py 


Jea = —km VP; + RT 


y= Jew (74) 


Applying this model to the wood fiber using the constant 
properties (without influences of humidity and temperature), 
Berger et al. [73] used a program translated with MatlabTM 
to obtain the curves that defined well the experimental data 
presented. Based on Whitaker’s model, Perré and Turner [83] 
presented a model applied to the wooden wall byRafidiarison 
[84]. The medium is supposed to be partially saturated, and 
mass transfers are applied in each phase. This model is given 
by Simo-Tagne [78]. 
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a > * * (75) 
22b +0 (pus) =m+mp 
ð = > * 
2 iT p p) = —Mp 
Py Uy = Py Ug — PgDer sV (2) 
PpUpy = —pcDpv (2) 
> Kyo ú (76) 
Ug = — 7 V(po) 
g 
> = _kig 
u i Vp) 
Pi = Pg — Pc (77) 


ply © + Ah, (m+ my) + hemp — poup TCh) + [(puu + 
Poup) Cpi + Zizan (pf UiCpi)|-VT = V (Aer VT) (78) 


Using experimental data given by Rafidiarison [84] and 
Rafidiarison et al. [85], the literature shows that this model 
gives satisfaction to predicting relative humidity, 
temperature, and moisture content of the wooden building 
walls in a temperate climate. 


5. Discussions 

The comprehensive bibliography of old and more recent 
hygrothermal transfer models for various building walls is 
reviewed in this paper. The hygrothermal transfer is very 
important for the design of a building envelope for thermal 
comfort and economic and energy analysis of the building. 
Several numbers of energy and mass conservation equations 
with different boundary conditions and input considerations 
have been presented in this paper for soil-based and wooden 
building walls. Some of the models are easier to use than 
others, though this doesn’t make them give better results 
because, in most cases, the validation with the experimental 
results has more errors. For example, the model of Lukoiv 
poses a lot of challenges because the conversion factor of 
water from a liquid state to a vapor state is located between 0 
and 1 due to the fact, they require the conversion factor to be 
obtained through a cumbersome process of the inverse 
method. Again, Duforestel model is easier to use than the one 
of Philip and De-Vries and can be applied to non-hygroscopic 
material in building, but the fact it did not take into account 
the gradient of moisture content is often the reason for the 
difference between experimental data and numerical data. 
The accuracy or otherwise of this model depends on the 
establishment of the right boundary conditions. Most of the 
research ignored the effect of hysteresis in their models, while 
very few considered the flow pattern of fluid through the wall 
surfaces. A literature review shows that most models were 
linked to Luikov equations for heat and mass transfer, and 
solutions were mostly by finite difference methods and 
implicit scheme of Crank-Nicolson to generate the numerical 
results. The thermophysical parameters of the medium were 
always considered constant during the process, but it will be 
good to consider swelling and shrinkage over time in the 
model as the material is influenced by environmental 
conditions. For example, in a high moisture environment, 
certain wood absorbs moisture and expands if used as a wall, 
while some shrink during the winter period. According to 
Ndukwu et al. [86], materials like wood, concrete or bricks 


40 


MC Ndukwu et al. /Future Technology 


used as a base material in wall construction if subjected to 
outdoor condition, undergoes continuous drying under 
ambient conditions. These materials, because they are 
hygroscopic, desorb or absorb moisture unless an 
equilibrium condition is maintained between the indoor and 
outdoor air, or the wall material surface is impermeable [87]. 
The subjection of these materials to external heat load from 
solar radiation and other high-temperature-producing 
sources like industrial activities results in continued moisture 
modification. Thus, moisture desorption or adsorption of 
most materials is a continuous process for the lifetime of the 
walls. The repercussion is the continuous modification of the 
thermophysical properties of this material. These require 
consideration in modeling the hygrothermal transfer for 
building walls. Therefore, as suggested by Ndukwu et al. [86], 
long-term behavioral models that will incorporate the 
continuous moisture loss or gain from building materials 
during structural application require investigation. The 
model should consider the effect of seasonal variations 
(winter, spring, summer, and autumn) and long-term 
meteorological data. Validations of the models showed the 
influence of wall, thickness, the density of the material, and 
climatic variations on the temperature and moisture 
evolutions within the building materials. Imaging models 
using software like COMSOL multi-physics, CFD etc. are scarce 
in the review bearing in mind that microscopic imagery is 
now deployed to measure the heat and moisture evolution in 
materials. Future models should include shrinkage or 
expansion influence, especially in fibrous materials like wood, 
as they respond to ambient variations. 


6. Conclusion 

The applications of various materials in building 
structures have been studied extensively. The study presents 
different models applied for predicting hygrothermal transfer 
for various building walls. Energy and mass conservation 
equations were applied with different boundary conditions, 
and thermophysical properties were presented in this paper 
for concrete, bricks, and wooden walls. Luikov models formed 
the basis of most models for porous materials. The 
parameters considered in most models were the nature of the 
materials of the wall, building orientation, variation in 
climate, the thickness of the wall, temperature, moisture 
changes, and the density of the material. Literature 
presenting imaging models using imagery software like 
COMSOL multi-physics, CFD etc. were scarce, considering that 
microscopic imagery is now deployed to measure the heat 
and moisture evolution in materials. Future models should 
include shrinkage or expansion influence on the fibrous 
material like wood due to their behavior under 
environmental conditions. 
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